Disentangling single-particle gap by electronic Raman absorption in electron-doped 

cuprates 
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In the under- to optimal-doping regimes of electron-doped cuprates, it was theoretically suspected 
that there is a coexistence of superconducting (SC) and antiferromagnetic (AFM) orders. The 
quasi-particle excitations could be gapped by both orders, and the effective gap is non-monotonic 
d-wave-like in the momentum space. Alternatively the gap was also speculated as a pure pairing 
gap. Using an effective microscopic model, we consider the manifestation of the quasi-particle gap 
in the electronic Raman spectra in a range of doping levels, where the relative strength of the SC 
and AFM order parameters varies. We demonstrate that from the electronic Raman spectra the 
effective single-particle gap can be disentangled into contributions from the two distinctive orders. 
This would help to tell whether the non-monotonic gap is due to the coexistence of SC and AFM 
orders. 



Understanding the pairing symmetry poses a strong 
constrain on the underlying superconducting mechanism. 
In hole-doped cuprates a consensus is reached that the 
predominant pairing channel is d x i_ y i [1]. The situ- 
ation is less clear in electron-doped cuprates, such as 
Nd 2 _ :E Ce x Cu04 and Pr 2 _ K Ce ;r Cu04. Tunneling [2] and 
early Raman measurements [3] suggest an s-wave pair- 
ing order parameter, while phase-sensitive measurements 
consistently imply <i-wave pairing[4], the latter of which is 
also supported by angle-resolved photo-emission spectra 
(ARPES) [5] and some penetration-depth measurements 
[6]. It was also speculated that there may be a transi- 
tion from d- to s-wave pairing as electron doping level 
increases [7] or as temperature dccrcases[8]. Even in the 
e?-wave picture, the situation is more complicated than 
that in hole-doped case. The ARPES [9] and the Raman 
scattering [10] experiments on Pro.sgLaCeo.nCuC^ and 
Ndi.85Ceo.i5Cu04, both of which are at the optimal dop- 
ing, revealed a non-monotonic ci :z ,2_ y 2-wave quasi-particle 
gap, i.e., the maximal gap value occurs between the nodal 
and antinodal directions instead of at the antinodal di- 
rection. On the other hand, a nontrivial evolution of 
the Fermi surface (FS) was found in electron-doped sam- 
ples [11, 12]. An electron-like FS pockets exist near the 
(±7r/4, ±7r) and (+7T, ±7r/4) regions in the momentum 
space for all doping levels in the range 0.10 < x < 0.15, 
while a hole-like FS pocket around (±7r/2, ±7t/2) does 
not emerge untill x = 0.13. Accordingly, Raman scatter- 
ing experiments revealed that the relative position of the 
B\g and Bi g peaks changes with doping. [13] In the opti- 
mally doped region, the Bi g peak appears at a higher fre- 
quency than the B\ g one, while in the overdoped region, 
it appears at a lower frequency than the B\ g peak. This 
was argued as due to the change of the shape of the single- 
particle gap that originates purely from superconducting 
pairing [10]. This picture is however not sufficient to ac- 
count for the evolution of the FS [11, 12], and is also in 
contradiction to the tunnelling measurements that con- 
sistently report an s-wave feature [2] . Another difficulty, 
which is perhaps common in both hole- and electron- 



doped cuprates, is the large signal in the A\ g channel of 
the Raman spectra [14], which would have been strongly 
suppressed by the screening effect in any one-band the- 
ory. A phenomenological two-band theory with tunable 
independent pairing gaps on the two bands seems to ac- 
count nicely the experimental feature [15]. Apart from 
the A\ g puzzle, the microscopic origin of such a theory 
is yet lacking. Alternatively, the non-monotonic gap was 
speculated as arising from the coexistence of d-wave su- 
perconducting (SC) and antiferromagnetic (AFM) orders 
[16, 17], even though the SC gap itself is in a typical 
monotonic d-wave form. In the under-doped regime, the 
hole-pockets do not appear yet, and the single particle ex- 
citations would be gapped everywhere in the momentum 
space. This is then not withstanding the s-wave feature 
in tunnelling[2] and specific heat measurements [8]. An 
issue does arise, however. ARPES and tunnelling only 
probe single particle excitations, it is not clear how the 
nontrivial gap manifests in two-particle properties and 
whether one can disentangle the single particle gap into 
contributions from different origins. This is the topic of 
the present paper. 

We calculate the electronic Raman spectra in the 
SC+AFM picture. In a range of doping levels, the rel- 
ative strength of AFM and SC orders varies, and the 
Raman spectra reveal distinctive peaks that can be as- 
sociated to the SC and AFM orders. We also compare 
the normal state with AFM alone. It turns out that 
the B2 g Raman channel is only sensitive to the SC gap, 
while the B\ g channel responses to both SC and AFM 
orders. In the presence of a hole-pocket, the B 2g channel 
reveals further a double peak structure, reminiscent of 
the two sets of normal state FS. Moreover, in all cases 
that we consider, the SC-related Raman peak in the B 2g 
channel is at lower energy than that in the B\ g channel, 
although the offset becomes much smaller in the coex- 
isting SC+AFM state. As an alternative to elastic neu- 
tron scattering measurements (not available yet), these 
results could help identify the mechanism of the non- 
monotonic gap and also whether SC+AFM coexists in 



2 



electron-doped cuprates. 

As a handy effective model, we consider the two- 
dimensional t-t'-t"-J model hamiltonian, 

H = -t ( C *W + h-c) - t' feW + h.c.) 

-t" E (4^+t.c.) + J^Sr4 (1) 

where (27)2, denote the nearest, second- 

nearest, and third-nearest neighbor bonds. For electron- 
doped cuprates, this model is defined in the hole picture 
so that c£. creates a cr-hole, and no double hole occupancy 
is implicitly assumed. A hole occupancy of 1 — x corre- 
sponds to a physical electron doping level x. Throughout 
this work, we use the magnitude of the nearest neighbor 
hopping integral as the unit of energy, so that t = — 1, 
t' = 0.32, t" = -0.16, and J = 0.3. As usual, the no 
double occupancy of holes is treated at the slave-boson 
mean field level where c a is replaced by \fxf a in the hop- 
ping terms, where f a is the fermionic spinon annihilation 
operator, and the spin exchange term is decoupled in a 
standard fashion as S t ■ Sj -> -(3/8)[x* ; E CT fLfj° + 

A^Cfa/it ~ fofil) + h.c] + (1/2)K £ CT trfUu, + * - 
j) + (3/8)(|xij| 2 + |Ay| 2 ) — miinj, where consistency re- 
quires Xij = T,Mlfjv)' A v = (fnfn - fnfji), and 

m i — J2a a (fLfic)/2- A different setting was used for 
the decoupling of the spin exchange term in ref. [16], al- 
though one does not expect qualitative changes in the 
results. Anticipating the uniform AFM order and d-wave 
X, Ajj = ±A for x- and y-bonds 



particle Hamiltonian, 



( A 1 e k 

£ k A 2 A k 

A k A 3 

\A k -e k 



A k 




(3) 



p + 2Jm, 



where A\ = e' k — p — 2Jm, A 2 = e 
A 3 = -e' u + fj, - 2Jm, and A 4 = -£ k + p + 2 3m. The 
single particle hamiltonian hk can be easily diagonolized, 
yielding a set of four eigenvalues i?k,n (n= 1, • • • , 4) and 
the corresponding eigenvectors |k, n) for each k. The 
single-particle 4x4 matrix Matsubara greens function 
G(k, iu) n ) = l/(iu; n l — /i k ), where 1 is the 4x4 identity 
matrix, can be expressed in terms of the eigen states as 



G(k,iiv n ) = E 



|k,n)(k,n| 

iul n — i?k,n 



(4) 



which forms the basis for further calculations. 

The electron Raman scattering measures the spectral 
function of the fluctuations of the effective charge density 
p = X^ kcr 7kC kcr Ck<T where is the Raman vertex that de- 
scribes the second order coupling between electrons and 
photons [18]. In our case the bare Matsubara propagator 
for p can be written as, 

X 11 {iv n ) = -— E Tr[7kG(k,iw„)7 k G(k, iw n + w n )].{b) 

Here T is the temperature, N is the number of lattice 
sites, and 



SC order, we set Xij 
respectively, and rrii — ±m for A and B lattice sites re- 
spectively. It is convenient to redefine a a = f a on the 
A-sublattice and b a = f a on the B-sublatticc. The mean 
field Hamiltonian of Eq.(l) can then be written as, in the 
momentum space, 

H M F = E[ £k ( a k<r &k <7 + nX -) + ( £ k - M)(aLr a k<r + & LA<t 
k,<r 

-2Jm ( r(aj C(T ak <T - & ko .&k CT )] 

-E A k( a kT&-ki +fok|a-ki +h.c), (2) 



7 k = (rii • V k )(n s • V k )/i k |j_ 



(6) 



is the matrix form of the Raman vertex, where rij and 
n s are unit vectors for the polarizations of the incident 
and scattered lights. Note that the spins do not couple 
to light to the leading order, so that the spin exchange 
term does not contribute to the Raman vertex, which is 
reflected in the formal replacement J — > in the above 
definition of 7 k . (We do not consider the higher order 
two-magnon Raman absorptions here.) The summation 
over iuo n in Eq.(5) can be performed analytically, and we 
obtain 



where e k = (—2tx — 3Jx/4)(cosfc x + cosfc y ), e k = 
— At'x cos k x cos k y — 2t"x(cos2k x + cos2fc y ), A k = 
— (3 J/4)A(cos k x — cosfcj,). Here p is the chemical po- 
tential that fixes the occupancy and the wave vector 
k is restricted in the magnetic Brillouine zone (MBZ). 
The order parameters are calculated self-consistcntly 
at each doping levels. For better convenience, we 
rewrite the mean field hamiltonian compactly as Hmf — 
SkV'k^kV'k, where Vk = (<XkT,&kT, a -ki> & -kj.) T are 
Nambu- Anderson four-spinors, and /i k is a 4 x 4 single- 



f(E m ) - f(E n ) 
L iu n + E m — E n 



k.rn.n 

xTr( 7k |k, m)(k, m| 7k |k, n)(k, n\)], (7) 

where / is Fermi-Dirac distribution function. The Ra- 
man spectral function is then given by 



R(E) = ImxOn 

7T 



E + i0+), 



(8) 



where E is the Raman shift. Note that theoretically 
R{E) is odd in E and must vanish at E = 0. For the 
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Big and B 2g channels, the Raman vertices are given by 

7k 19 = ("fk X ~ 7k S ')/2 an( i 7k 2 ' = 7k^' respectively. The 
B 2g channel is d^-symmetric and mainly probes quasi- 
particle excitations in the nodal directions, while B\ g is 
d x 2_ y 2 symmetric and mainly probes excitation in the 
antinodal directions. These properties enable Raman 
scattering to selectively probe excitations in the momen- 
tum space, and is applied in hole-doped cuprates to verify 
the d-wave nature of the pairing gap. [18] Both B ig and 
B 2g vertices are odd in parity, the Raman absorption in 
these channels are not re-normalized by Coulomb inter- 
actions. In contrast, in the fully symmetric A\ g channel, 
the re-normalization is severe and depends on the de- 
tailed band structure (or the Fermi surface harmonics). 
Due to such complications, we shall concentrate on the 
simpler B\ g and B 2g channels in this work. 

Before the exposition of Raman spectra in the 
AFM+SC states, it is instructive to understand, for com- 
parison, the Raman spectra in the non-superconducting 
states. First, in the normal metallic state with neither 
AFM nor SC, no Raman spectra are anticipated un- 
less scattering sources are introduced. This is because 
no finite-energy particle-hole excitations with zero mo- 
mentum transfer can be achieved in ideal metals. Sec- 
ond, in the normal state with AFM, the band is split 
into upper and lower bands. Particle-hole excitations 
vertically across the sub-bands is possible, and one ex- 
pects nonzero Raman absorption. In this case we show 
that the Raman response in the B lg channel scales 
with the square of the AFM order parameter, while the 
B 2g channel is completely insensitive to AFM for the 
ranges of doping under concern. In the AF state, the 
Hamiltonian in momentum space can be simplified as 
Hmf = Ska^L&^ka, where 7y k(T = (a kcr , 6 kcr ) T , and 
Ck.cr = (e k — aO^o + £kCi — 2Jmaa 3 , where <jq is the 
unit 2x2 matrix and <j\^ are Pauli matrixes. The spin- 
dependent single particle Green's function is given by 
<? a (k, iuj n ) = l/(iu> n ao — £k,o-), which we rewrite for later 
convenience as, 



g a (k, iuj n ) 



2 ^ i. 



jt±l + M — £ k — vE k 

e k cri - 2Jmaa 3 
x K + v — ), 



(9) 



where Ek = yej + 4J 2 m 2 . The Raman vertex is given 
by 7k = (n» • V k )(n s • V k )£ ki(T )l J^o- By explicit algebra, 
we obtain 7^ 2s = a k £r 0j where a k = — At'x sin k x sin k y , 
and 7]f l9 = & k oo + c k <ri, where & k = At"x(cos2k x — 
cos2k y ) and c k = tx(cosk x — cosk y ). In each spin- 
channel, the Raman response is given by a formula simi- 
lar to Eq.(5). The total response function summed over 
spin species is then given by 

f(u 2 E^ + e k - n) - /Oi-E k + e k - y) 



X 



*9 = J_ V" 

4A ^ 
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FIG. 1: Quasiparticle gap A 9 and Raman spectra at the dop- 
ing level x — 0.143. (a) A 9 as a function of 9 (in degrees) 
in the AFM normal state, (b) Raman spectra corresponding 
to (a), (c) A 9 as a function of 9 in the AFM+SC state, (d) 
Raman spectra corresponding to (c). 



£ k cri - 2Jmaa 3 
xTr[a k (cr + v\ = ) 

e k <7i - 2Jmoo 3 
xa k (cr + v 2 J, 



(10) 



where V\,v 2 = ±1. The trace term gives 2a 2 . (1 + V\V 2 ), 
but only the cases with V\ = — v 2 contributes eventually 
due to the two fermi functions in the first line. Therefore, 
\ B2g = in the normal state even with AFM order. On 
the other hand, for the B lg channel, the final result is 

o 72 2 

— jy- + e k - M ) - f(E k + e k - M )] 

1 1 



,Bi 



x-tt( — 
El y iv n 



2E V iv n - 2Bh 



-), 



(11) 



k,cr,^i,z/2 



which scales with m 2 as we emphasized. We emphasize, 
however, that the conclusion that B 2g = depends on 
the bare band structure we use. In principle, including 
4-th and further neighbor hopping would cause a nonzero 
B 2g , but we expect it to be much weaker as compared to 
Bi g on general grounds. 

We now present and discuss the numerical results. The 
case of x = 0.143 is presented in Figs.l. In order to 
understand the Raman spectra we present the quasi- 
particle gap A g as a function of the angle 9. Henceforth 
9 is defined as the angle between k — Q and — Q where 
Q = (it, 7r). The gap A g is then defined as the minimum 
of the quasi-particle energy in the momentum space along 
a cut of k at the same angle 9. Clearly A s (0) measures 
the gap in the nodal direction, while A g (45°) measures 
that in the antinodal region, given the fact that the nor- 
mal state fermi surface is close to the antinodal point 
(•7T, 0). We begin with the AFM normal state (by sim- 
ply erasing the SC order from a state with AFM+SC). 
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The gap is presented in Fig. 1(a). For < 9 < 10° and 
25° < 9 < 45°, the normal state gap A g = 0, roughly re- 
flecting where fermi surfaces appear. For 10° < 9 < 25°, 
the normal state is gapped by the AFM order, with a 
non-monotonic gap A g as a function of 9 due to the un- 
derlying band structure. The Raman spectra in this state 
are presented in Fig. 1(b). As we analytically proved, 
for the model at hand, there is no B 2g Raman absorp- 
tion, but the Big absorption is strong, and develops a 
threshold Raman shift roughly twice of the maximal gap 
in Fig. 1(a). Now turning on the SC order, the normal 
state gap superimposed by the monotonic d-wave pair- 
ing gap leads to the gap structure presented in Fig. 1(c). 
Similar result was also reported elsewhere. [16] The max- 
imum of the gap occurs at = 18°. In the experimental 
case[9, 10], 9 = 20° and 9 = 15° in Pr .89LaCco.iiCu0 4 
and Ndi.85Ceo.i5Cu04, respectively. The corresponding 
Raman spectra are presented in Fig. 1(d). In the B lg 
channel, the higher energy peak can be associated to 
the AFM normal state Raman peak but slightly shifted 
due to pairing. The new lower energy peak is definitely 
caused by the SC order, as seen from the fact that the 
energy of this peak is twice of the pairing gap at the 
anti-nodal direction. The B 2g channel, which is silent in 
the normal state, becomes active in the AFM+SC state, 
and the peak energy is therefore entirely determined by 
pairing alone. Interestingly, because of the two kinds of 
fermi surfaces in the case under concern, two absorption 
peaks appear in the B 2g channel. The energies of the 
peaks are roughly (but smaller than) twice of the val- 
ues of the pairing gaps at the two fermi surfaces. This 
is consistent with the general observation that the B\ g 
Raman vertex is maximal in the the anti-nodal region, 
whereas B 2g vertex is maximal in the nodal region and 
therefore does not see the full pairing gap. Summarizing, 
we claim that electronic Raman scattering not only dis- 
entangles the different contributions to the total single 
particle gap, but also tells the number of distinct fermi 
surfaces, given the monotonic d-wave pairing. Finally, 
the SC-driven B 2g Raman peaks are at smaller energy 
than the SC-driven B\ g peak energy, even though the 
single particle gap is as non-monotonic as in Fig. 1(c). 
In experiment [10], there appears to be only one peak in 
both Big and B 2g channel, and the B 2g channel is situ- 
ated at a larger energy. It is also not clear whether an 
AFM-driven Raman peak is present experimentally. 

We now shift to the under-doped AFM+SC regime 
where AFM is even stronger. For x = 0.13, the AFM 
normal state single particle gap is presented in Fig. 2 (a). 
Now the gap vanishes only near the anti-nodal region, re- 
flecting the fact that the hole fermi pocket is absent and 
only one kind of fermi surface appear. The correspond- 
ing Raman spectra are presented in Fig. 2(b), where again 
the B 2g channel is silent and the B\ g channel develops 
AFM-driven peaks. The additional peak at higher energy 
arise from the van-Hove singularity in the band structure. 
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FIG. 2: The same plot as Fig.l but for x = 0.13. 
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FIG. 3: (a) A 9 as a function of 9 in the SC state at the doping 
level x = 0.15. (b) Raman spectra corresponding to (a). 



Turning on the SC order, the single particle gap is pre- 
sented in Fig. 2(c). Due to the underlying AFM normal 
state gap, the effective gap here is already nonzero in the 
nodal direction, and is overall non-monotonic. The cor- 
responding Raman spectra are presented in Fig. 2(d). In 
the Bi g channel an SC-driven new Raman peak appears 
in addition to the AFM-driven peaks at much higher en- 
ergy. The new peak is at an energy twice of A g = 0.02|t| 
at 9 — 45°. On the other hand, the B 2g channel shows a 
Raman peak with Raman shift very close to but slightly 
smaller than the energy of the SC-driven peak in the Big 
channel. Since the B 2g channel is silent in the normal 
state, we conclude that the gap seen by the B 2g channel 
has nothing to do with the nodal effective gap, but is 
rather determined by the pairing gap. The single-peak 
structure seen in the B 2g channel is also consistent with 
the the absence of hole pockets. A disentanglement of 
single particle gap is therefore again possible. 

In the over-doped region, for example at x = 0.15 in 
our case, the self-consistent calculation yields that AFM 
order disappears, and we recover a single-band with SC. 
The situation is then formally similar to the hole-doped 
case. For completeness we present the quasi-particle ex- 
citation gap Ag and the corresponding Raman spectra 
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in Fig. 4. The gap as a function of 9 is typical of a pure 
d-wave pairing gap. The B\ g channel develops an ab- 
sorption peak at exactly twice of the maximal pairing 
gap, while the ratio is roughly 1.5 in the B 2g channel. 
Because of the absence of AFM order, both channels 
develop single peaks. The relative peak position is in 
agreement with experimcnt[13], but the relative intensity 
of the peaks is not, possibly due to more fine details of 
the materials. For example, inter-band scattering from 
the oxygen p-band to the upper Hubbard band under 
concern may modify the effective mass of the electrons, 
leading to photon-energy dependent changes in the effec- 
tive Raman vertices. This effect has been experimentally 
observed. [10] 

In conclusion, we worked out the consequence of AFM 
order coexisting with the SC order in the single-particle 
gap as well as the electronic Raman scattering. In par- 
ticular we show that using electronic Raman scattering it 
is possible to disentangle the effective single particle gap 
into distinctive contributions from AFM and SC orders. 
Combining the studies in the one-band picture (but with 
non-monotonic pairing gap) [16] and the phenomenolog- 
ical two-band picture (with independent pairing gaps on 
the two bands) [15], and by comparing to existing and 
forthcoming experiments, a strong constrain can be made 
on the AFM origin of the non-monotonic gap as well as 
the evolution of the fermi surface with doping in electron- 
doped cuprates. 
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